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Abstract 

CN ■ The spectrum of complex PT-symmetric potential, V(x) = igx, is known to 

CN ' be null. We enclose this potential in a hard-box: V(\x\ > 1) = oo and in a soft- 

> : 

Q\ , box: V(\x\ > 1) = 0. In the former case, we find real discrete spectrum and 

' the exceptional points of the potential. The asymptotic eigenvalues behave as 

On i „ 

E n ~ ti . The solvable purely imaginary PT-symmetric potentials vanishing 

asymptotically known so far do not have real discrete spectrum. Our solvable 



soft-box potential possesses two real negative discrete eigenvalues if \g\ < 
(1.22330447). The soft-box potential turns out to be a scattering potential 



| not possessing reflectionless states. 

CT PACS NO. : 03.65.Ge 
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A complex Hamiltonian that does not change under the joint transformation of Parity 
(P: x — > —x) and Time-reversal (i — > —i) is called PT-symmetric: (PT)H(PT) -1 =H or 
[PT,H]=0. When two (linear) operators A and B commute, we can arrage for simultaneous 
eigenstates. However, if one of the operators is anti-linear: A (cip) = c*ip*, then for A and 
B we will have simultaneous eigenstates if only the eigenvalues of B are real. In case of 
complex eigenvalues of B, the simultaneity of eigenstates does not hold. Further, it follows 
that B can have real or complex conjugate pairs of eigenvalues. Therefore, PT being anti 
linear a Hamiltonian commuting with it may have entire discrete spectrum as real or as 
complex conjugate pairs or both. In the first case PT-symmetry is unbroken or exact in the 
second case it breaks down spontaneously. 

The revelation that a complex PT-symmetric potential [1] too may possess real discrete 
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spectrum has lead to very interesting developments [2]. In this regard, the exactly solv- 
able potentials have been very revealing [3]. Quasi-exactly solvable complex PT-symmetric 
potentials [4] have also been helpful wherein one gets only a certain number of eigenvalues 
and eigenstates witb some restriction on the parameters of the potential. Next some sim- 
ple semi-analytically [5] solvable potentials have been proposed wherein one gets implicit 
equation for energy to calculate the discrete spectrum iteratively and numerically. These 
potentials have been useful and handy in some tedious calculations 

For the potential Vb(x) = — (ix) v [1], a few critical values of v in (1.42207,2) are known 
where the spectrum has a few real eigenvalues and for v > 2 the entire spectrum is real. 
For analytically solvable potential, Vs(x) = — Visech 2 a; + iV 2 secha;tanha;, [3] the discrete 
spectrum is real if | V2 j is less than or equal to a critical value yJVi + 1/4. Otherwise the 
spectrum consists of complex conjugate pairs of energy. The semi-analytically solvable 
potential Vz(|x| < 1) = iZ9(x),V z (\x\ > 1) = 00 [5] has been found to have critical 
values of Z above which PT-symmetry is spontaneously broken. Here 9(x < 0) = — 1 and 
6{x > 0) = 1. 

Nowadays, with the advent of the packages like Mathematica one can perform quick and 
almost error-free calculations involving even the higher order functions. This gives us a 
broader scope to extend at least the class of semi-analytically solvable potentials. In this 
Letter, we present the complex PT-symmetric potential V(x) = igx in two modifications. 
These two models are semi-analytically solvable which also arise due to the following open 
questions in PT-symmetric quantum mechanics [6,7]. 

Another way of visualizing real spectrum of a non-Hermitian Hamiltonian is through the 
exceptional points of a Hamiltonian as proposed by Kato (1960) [8,9]. Let the eigenvalues 
E n (\), of H\ — Hi + \H 2 be analytic function of A, where Hi, H 2 are Hermitian Hamilto- 
nians. If A is analytically continued say as A = iji in a complex domain, then at certain 
values of fj, — /i c two real distinct eigenvalues and eigenstates may coalesce (merge) in to 
one. Further, if /i > /i c the coalesced eigenvalues reappear as complex conjugate pair. In 
this complex extension of the parameter A the Hamiltonian becomes non-Hermitian. We 
strongly feel that PT-symmetry or more generally pseudo-Hermiticity [10] (rjHi]^ 1 = W) 
may be the necessary condition for the existence of exceptional points in R ifM . 

So far we do not know the necessary and sufficient condition to determine a priori whether 
a given complex PT-symmetric potential will have real discrete spectrum. However, case 
of V(x) = igx is well studied wherein using the arguments of Stokes and Anti-stokes lines 
the spectrum has been found to be null. An interesting a priori and necessary condition in 
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this regard is that the potential ought to have a pair of classical turning points (solutions of 
E — V(x)) of the type ±a + ib [11] at a given real energy. In addition to this, any singularity 
of the potential should not lie on the line from z = —a + ib to z = a + ib. This explains why 
the real discrete spectrum of Vl(x) = igx[l],VE(x) = iV 2 tanhxsechx, Vq{x) = iV 2 xe~ x2 
and Vc(x) = iV2x/(l + x 2 ) for any real value of V 2 is empty. Since only Ve(x) among 
these potentials is analytically solvable having the property of shape-invariance, therefore 
an argument based on supersymmetry may be ruled out for the emptiness of the real discrete 
spectrum. However, when a real attractive potential of sufficient strength is added, these 
potentials [3] do possess real discrete spectrum (e.g., Vho( x ) = x 2 + igx). 

In this work, we study the purely imaginary potential, Vl(x), within a hard-box ((V(|x| > 
1) = oo) or a soft-box (V(|x| > 1) = 0) for real discrete spectrum. The abovementioned 
potentials, Vg(x) and Ve(x) may also be seen as purely imaginary PT-symmetric potentials 
vanishing asymptotically and possessing the real discrete spectrum as null. Similarly, we 
find that simple purely imaginary potential V(x) = iV 9(x) placed within the soft box also 
has real discrete spectrum as null. This gives us an additional motivation to discuss Vl(x) 
placed within a soft box. 

First we study 

V(x) =igx,V(\x\ > 1) = oo. (1) 

Taking 2m = 1 = h 2 in one-dimensional time-independent Schrodinger equation, we get the 
wave solution as tp(x) = AAi(z) + BBi(z) with z = (E — igx)/q, q = {g 2 ) 1 ^. The Dirichlet 
boundary conditions that ip(x = ±1) = and their simultaneous consistency gives rise to 
an implicit and transcendental equation for E as 

Ai[(E + tg)/q]Bt[(E - tg)/q] - Ai[(E - tg)/q}Bt[(E + tg)/q] = (2) 

Here Ai(z) and Bi(z) are the Airy functions. The eigenfunction is expressible as 

i/>(x)=C(E,g)Ai[(E + igx)/q] + Bi[(E + igx)/q], -1<x<1, (3) 

where 

Bi[(E + ig)/q] _ Bi[(E-ig)/q] 

Ai[(E + ig)/q] 1 > 9) ~ Ai[{E-ig)/qy 1 ] 

For the real branch of q, the Eq. (2) say f(E) = is a real equation of E on real line. 
Therefore as per the fundamental theorem of algebra the roots (eigenvalues) will be either 
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purely real or complex-conjugate pairs. In the former case, ip(x) will also be an eigenstate 
of the anti-linear operator PT (x — > —x and i — > — i). In the latter case, ip(x) will no 
more be the eigenstates of PT and the PT-symmetry will be spontaneously broken. We 
employ Mathematica to perform the calculations. We find the first critical value g± = 
12.3124556046 of g. If \g\ < gi, each calculated eigenvalue (root of from Eq. (2)) is real 
and discrete. Then if \g\ = g 1: we find that only the lowest pair of real discrete eigenvalues 
merge (coalesce) into one real eigenvalue to make way (when \g\ > g\) for a single pair of 
complex conjugate eigenvalues signifying spontaneous breakdown of PT-symmetry. When 
\g\ = gi, we notice the merger of lowest real eigenvalues at E — 7.1086. A sudden jump 
discontinuity is observed in the ground level while crossing a critical value e.g., Eq{Y2.?>1) = 
7.0316, f? G (12.32) = 21.7209. The first five eigenvalues for g = 12.31 are 7.03165, 7.1848, 
21.7217, 39.1884, 61.4929. On the other hand, the first five eigenvalues for g = 12.32 
are 7.1097 ± .1342i, 21.7209, 39.1880, 61.4926(61.6850). The number in bracket equals 
257r 2 /4. The next critical values of g are g 2 ~ 53.18, g 3 xs 122.90, ... The first five eigenvalues 
for g = 53 are 16.4942 ± 24.431H, 29.7350, 31.6153, 58.2078 and for g = 54 they are 
16.7009 ± 25.0725i, 30.8513 ± 1.973H and 58.0801. 

Thus each critical value of g indicates the onset of the removal of lowest pair of real 
discrete eigenvalues. Removal of only ground level of a potential from its supersymmetric 
partner potential has earlier been observed in supersymmetric quantum mechanics [12]. In 
contrast to this, a complex PT-symmetric potential having infinite real discrete spectrum 
manifests in the removal of the lowest pair of real eigenvalues as its parameter passes through 
several critical values. When this happens the ground state eigenvalue becomes discontin- 
uous and we have A(gi) = [(-^(M = g\ + e) — E G {\g\ = g 1 — e)]. When the parameter g 
increases the net density of states of real eigenvalues reduces. 

Seeing the potential (1) as a non-Hermitian PT-symmetric perturbation to the one di- 
mensional box having E® = n 2 ir 2 /4 presents interesting features. We find that there are two 
remarkable features of the complex PT-symmetric perturbation (1) in contrast to its real 
Hermit ian counterpart when V(|x|<l) = ^x. In the former the value of g increases 

the lowest eigenvalues(approximately corresponding to E®) of the real discrete spectrum get 
removed reducing the number of real eigenvalues below a fixed energy E = E*. On the other 
hand, in the Hermitian case, these lowest eigenvalues get perturbed and get pushed down 
to become negative however the number of levels below E = E* does not change. Secondly, 
we find that the higher real eigenvalues behave as E n ~ E® + e\ whereas in the Hermitian 
case they are like E n = E® — e|, ei, e 2 are small real numbers. 
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We also find that there are bands of values of g, wherein for a fixed value of g, we 
get a continuous band of energies which satisfy the eigenvalue equation (2), however the 
corresponding eigenstates are null. The discrete eigenvalues do exist on the left and right 
of this energy-band. For example when g = 3.4, we get a continuous band of energies 
23.2 < E < 26.5 wherein every energy is an "eigenvalue" with a null eigenvector. For 
g = 5.0z, we get such an energy band as 32.8 < E < 34.2. 

The discrete spectrum of the Hamiltonian with a hard-box potential can also be con- 
veniently calculated by diagonalizing the matrix H rs =< r\H\s >, where \r > are 
the complete orthnormal eigenstates of the one-dimensional hard-box potential namely 
Vff(M < 1) = 0, > 1) = oo. We get for our hard-box potential 

r 2 TC 2 r , _., /9 lQiqrs r l + ( — 1)* „ „ , r 

^ = — ^ + (-l) f/2 ^ 2(r2 ^ 2)2 [ ^-], t = r + s + l,r,s = l,N. (5) 

Here N denotes the truncation or the order of the "infinite" dimensional matrix H r s . The 
eigenvalue determinant det \H — EI\ NxN is a very fast converging function of N. Conse- 
quently, the critical values of g discussed above hardly depend upon N. This in turn also 
means that first N eigenvalues (real or complex) calculated by fixing a value of g and the 
dimensions of the H matrix N x N would be quite accurate. We would like to remark that 
the method of diagonalization of eigenvalues supplements the calculations using Eq. (2) 
specially when eigenvalues are complex conjugate pair. 

The PT-symmetry is also conceived as P-pseudo-Hermiticity [9], it is therefore worth 
pointing that the matrix H r s (5) for the PT-symmetric matrix is pseudo-Hermitian for any 
value of g with parity operator given by P = (— l) r 5 rjS . Other property of H r s is that it 
is complex symmetric matrix. The experience gained in the diagonalization of a hard-box 
Hamiltonian is markedly different from that of diagonalization of the well known Hamiltonian 
of D. Bessis H DB = p 2 — igx 3 [1]. 

The diagonalization of Hob in the real harmonic oscillator basis by truncating the di- 
mensions of the matrix to finite: N x N, contrary to yielding the whole discrete spectrum as 
real for any value of g [1], yields [13] a finite number (TV) of eigenvalues as real if \g\ > g c (N). 
Here g c (N) is crucially decreasing function of N. So eventually as N — > oo, g c (N) vanishes 
justifying entire spectrum to be real for any real value of g. 

The complex PT-symmetric Potential 

V R (\x\ < 1) = -V, + iV 2 , V R (\x\ > 1) = (6) 
possesses scattering states for E > and bound states for E < 0. The interesting scattering 
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properties of (6) have earlier been studied for this potential [14]. We further find that when 
V\ = 0, the purely imaginary potential (6) has the real discrete spectrum as empty. This 
potential along with Ve(x) and Vg(x) [3] constitute purely imaginary and asymptotically 
vanishing PT-symmetric potentials that do not possess any real discrete eigenvalue. We now 
study the second modification of igx-potential 

V(\x\<l)= i9 x,V(\x\>l) = 0, (7) 

the soft-box potential that is exceptional in this regard. Assuming 2m = 1 = h 2 , we consider 
the Schrodinger equation for (7) when E < 0. Using the Airy functions again, we write the 
wave solution of Schrodinger equation in three regions as ip(x < —1) = ce KX , 1 < x < 
1) = aAi(z) + bBi(z),tp(x > 1) = de~ KX . Matching the wave solution at x — ±1 and 
eliminating a, b, c, d, we get 

[nqAi(s) + igAi'(s)] [nqBi(t) - igBi'{t)\ = [nqBi{s) + igBi'(s)] [KqAi{t) - igAi'{t)\ (8) 

Here s,t = (E ± ig)/q,z — (E — igx)/q,q = (g 2 ) 1 ^ and k = \J—E. This equation is an 
implicit equation of —E which is real for real negative values of E. So as per the fundamental 
theorem of Algebra it will have either negative real roots or complex conjugate roots. By 
the evaluation of the roots of Eq. (8), we find that if \g \ < 1.223830447 = g c there exist 
two negative real discrete eigenvalues. Further if \g\ = g c , we find that two real discrete 
eigenvalues merge at E = .24994. For the values of \g\ > g2, the eigenvalues are complex 
conjugate pairs. As \g\ is decreased one level is pushed down to negative values, the other one 
shifts closer to zero. For instance for \g\ = 1.2, E x = —0.40891, E 2 = —.14426. For \g\ = 0.6, 
we get Ei = -2.58012, £ 2 = -0.00275. When \g\ = 0.1, we get E x = -9.29466, £ 2 = 
— 1.7849 x 10~ 6 . We conjecture that the PT-symmetric potentials which have both bound 
and scattering states will have only one critical value of parameter as found here for (7) and 
in the models in Refs. [3]. 

The limit when g approaches zero is not transparent here. For this we use the asymptotics 
or the Airy functions Ai(s), Bi(s) ~ s3e T2s ^ /2 / 3 for large values of s xs t(~ ^573)- We find 
that the Eq. (8) becomes independent of energy and it is trivially satisfied for g = 0. This 
confirms the Hermitian free particle limit of the soft-box potential wherein any energy is a 
solution of Eq. (8) and hence discrete energies do not exist. 

Thus, we have at least one purely imaginary PT-symmetric potential vanishing asymp- 
totically and possessing real discrete spectrum. This investigation receives significance due 
to interesting results [15] of Hermitian quantum mechanics like: a real potential well such 
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that V(x)dx < oo has at least one discrete eigen value irrespective of the depth and 
width of the potential. In the absence of potentials like (7), it would be as though a purely- 
imaginary PT-symmetric potential that vanishes asymptotically does not possess real dis- 
crete spectrum. 

Earlier we have proposed generation of discrete eigenvalues of a Hermitian potential 
V(x) = —x A [16] by imposing PT-symmetric boundary condition on the wave function. At 
these discrete energies the reflection probability vanishes and potential becomes reflectionless 
[16,17]. If we analytically continue k as k — ik in the eigenvalue equation (8) we get 

[kqAi(s) + gAi'(s)][kqBi{t) - gBi'(t)} = [kqBi(s) + gBi'(s)][kqAi(t) - gAi'(t)], k = Ve, E > 0. 

(9) 

This equation is complex on the real line of E, so according to fundamental theorem of 
algebra it will have complex roots and some of these roots may or may not be real. For 
E > 0, the states are scattering states therein real discrete eigen value/ values if occur would 
be embedded in positive energy continuum. These states are discrete energy eigenstates 
generated by imposing PT-symmetric boundary conditions as ip(±x) ~ e ±ifcr , k = \fE ) E > 
instead of the usual ip(±x) ~ e ±KX , k = \/—E, E < 0. However, for our soft-box potential 
(7) we have carried out a very careful search and we do not find real solution(s) of the Eq. 
(9). Therefore, such reflectionless states do not seem to occur for this potential. 

One interesting common feature of gx-potentials lies in the presence of q = (— g 2 ) 1//3 (cube 
root of — g 2 ) which has one real (go) and two complex (gi,g2) branches. For the Hermitian 
case, we find that all the three branches of q give identical discrete spectrum. However for 
a given eigen value, we have three eigenstates ip qo (x),ip qi (x) and tp q2 (x). The first one is 
real and other two are complex valued functions of x. Since there can be no degeneracy in 
one-dimension, we therefore find that the last two eigenstates differ from the first one by 
multiplicative complex constants only. For the non- Hermitian case of igx-potentials (1,7), we 
find that the complex branches q± t2 may not yield the full real discrete spectrum as generated 
by the main (real) branch go- Only first few eigenvalues of the main branch are repeated. 
Amusingly, the continuous band of the higher energies turns out to be "eigenvalues" however 
with null eigenvectors. Every repeated real eigenvalue for the three branches once again 
yield three corresponding complex eigenfunctions: ijj qo (x),ijj qi (x) and ijj q2 (x). Only the first 
one is also eigenstate of PT, whereas other two are not. We once again find the last two 
eigenstates differ from the first one by a complex multiplicative factor. This wards off two 
fake occurrences: degeneracy (in one dimension) and the breaking PT-symmetry (despite 
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real eigenvalues) when one chooses a complex branch of q. The complex branches of q do 
not create any new energy eigenvalue other than those yielded by the real branch. 

Finally, we conclude that though the spectrum of V(x) = igx is null, however, the 
spectrum of two of its modifications are interesting. These are semi-analytically solvable 
models. The PT-symmetric hard-box potential helps studying issues like exceptional points, 
reality of eigen values, spontaneous PT-symmetry breaking and some interesting features 
of Hamiltonian-diagonalization. The reduction in the density of real eigenvalues as the 
parameter g increases is found to be disparate and strong feature of PT-symmetry. We 
conjecture that PT-symmetry or pseudo-Hermiticity could be a necessary condition for ob- 
serving exceptional points in the analytically continued complex Hamiltonian. Our soft-box 
potential is the first solvable model purely imaginary PT-symmetric which vanishes asymp- 
totically and has real discrete spectrum. In the light of this new possibility, it could further 
be interesting to search for purely imaginary and asymptotically vanishing potentials hav- 
ing discrete spectrum. It may be that the higher-order imaginary potentials of the type 
V n {x) = iV2(sech) n xtanhx, n > 1 do possess real discrete spectrum since they will have 
PT-symmetric turning points of the type ±a + ib Since these are not analytically solvable 
it would require numerical calculations. 

We do not find reflectionless states embedded in positive energy continuum in our soft- 
box potential. 

Very importantly, it turns out that all the three mathematical branches of q altogether 
output an unambiguous and non-anomalous physical result (spectrum) in both the cases 
whether we have Hermitian or PT-symmetric Hamiltonians. Perhaps only gx or igx poten- 
tials that are amenable to analytic solutions facilitate such a study. This study has not been 
without interesting surprises as mentioned earlier. Just recently both gx and igx potentials 
have been found very interesting for their complex classical trajectories [18]. 

We wish that the two modifications of igx-potential presented here will also be helpful 
in understanding other issues of PT-symmetric quantum mechanics in future. 
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